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We investigate the properties of the spectral function A(u, U) of correlated 
electrons within the Hubbard model and dynamical mean-field theory. Curves 
of A(u>, U) vs. <jj for different values of the interaction U are found to in- 
tersect near the band-edges of the non-interacting system. For a wide range 
of U the crossing points are located within a sharply confined region. The 
precise location of these "isosbestic points" depends on details of the non- 
interacting band structure. Isosbestic points of dynamic quantities therefore 
provide valuable insights into microscopic energy scales of correlated systems. 

PACS numbers: 71.27. +a 



1. INTRODUCTION 

A family of non-monotonic curves, obtained by plotting a quantity 
f(x,y) as a function of one of its variables (say, x) for different values of 
y, will in general intersect. The crossing points are located along a curve 
x*(y) denned by 

9f(x,y) 



dy 



= 0. (1) 

x*(y) 



In physics, chemistry and biology these crossing points are sometimes found 
to be confined to a remarkably narrow region, or even located at a single 
pointP thus leading to a conspicuous feature termed isosbestic point? 3 In 
the former case x* depends only weakly on y, while in the latter case x* does 
not depend on y at all. 

For example, the curves of the specific heat C(T,X) vs. temperature 
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T of numerous strongly correlated fermionic systems are known to cross 
once or twice when plotted for different values of a second thermodynamic 
variable xPln particular, crossing points are observed for different pressures 
(X = P) in normalfluid 3 He (Fig. ^i) 5 and heavy-fermion systems such as 
CeAi3 (Fig. ^)PBy changing the magnetic field (X = B) the same feature 
is seen in heavy-fermion compounds such as CeCvL^ x A\ x (Fig. [20 P 

Crossing points of specific heat curves are also observed in lattice models 
for correlated electrons such as the one-band Hubbard model 

H = J2 Uj c\ a c ja + n H n n -^^2n ia , (2) 

ija i ia 

where c| CT are creation operators for an electron at site i with spin a and = 
cj^Cjp are density operators. The model contains the hopping amplitude tij , 
the local Coulomb interaction U, and the chemical potential /i as parameters. 
At half-filling the curves C(T,U) vs. T always cross at two temperatures. 
This is observed, for example, in the case of nearest-neighbor hopping in 
d = lJHED d = 2 ED anc i d = oo (Fig. [IJi)Pl as well as for long-range hopping 
in d = lP 

As shown in Ref. [I] the existence of crossing points in the specific heat 
curves, and the fact that they may be quite sharp, can be linked (i) to a sum 
rule for the change of the entropy S(T, X) with respect to X in the limit of 
T — ► oo, and (ii) the properties of the susceptibilities x^(T, X) = d n £/dX n , 
where £(T, X) is the conjugate variable to X. Furthermore, the fact that 
the high-temperature crossing point in C(T,U) for the one-band Hubbard 
model occurs at a nearly universal value C*/ks — 0.34 was shown to be 
a consequence of the existence of two small parameters: the integral over 
the deviation of the density of states from a constant value, and the inverse 
dimension, 

The isosbestic points discussed above all appear in curves of the specific 
heat C(T, X) vs. temperature T when plotted for different values of another 
thermodynamic variable X. On the other hand, such points are known to 
occur also in dynamic quantities, e.g., in the optical conductivity a(u,n) 
of the high-T c material Nd2- x Ce a; Cu04_ 2/ where n is the density or doping 
(Fig. ^),^and in the Raman response xi - 1 ^) of the Hubbard modelPln 
general, there is no reason for the intersection of these curves to occur at one 
sharp frequency uj*, i.e., to be completely independent of any other param- 
eter defining the family of curves. There are, however, at least two classes 
of isosbestic points which are genuinely point-like. One is the exact crossing 
of curves described by a scaling function in critical phenomena; see, for ex- 
ample, the crossing of conductance curves near the Anderson transition!^ 
Another one is found in optical spectroscopic studiea^EI Q f systems con- 
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Fig. 1. Examples of crossing points in various correlated systems: (a) Specific 
heat of normal liquid 3 He measured at different pressures^ (after Ref. QJ. 
(b) Specific heat of CeAl3 measured at different pressures^ after Ref. HJ. (c) 
Specific heat of CeCus.sAuo.s measured at different magnetic fields^ (after 
Ref. QJ. (d) Specific heat of the Hubbard model calculated for different 
interactions^ (after Ref. H}. (e) Optical conductivity of Nd2- x Ce x Cu04_ y 
measured at different doping (after Ref . IT1J) . (f) Isosbestic point in ultra- 
violet spectra of three solutions of phenol and nitrate (after Ref. IT7)) . 
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sisting of two components with densities n a , rib an d n = n a + n& = const. 
In this case the absorbance as a function of frequency a(ui, n a ) depends only 
linearly on the density, i.e., has the special form 

a(u,n a ) = n a f a (oj) + (n - n a )f b {u). (3) 

If f a and f coincide at a frequency u>* , i.e., f a {u*) = fb(uj*), then 

da(uj,n a ) 
dn a 

This implies that for all densities n a the absorbance curves intersect at one 
frequency uj* (or the equivalent wave length, Fig. U)- Similar isosbestic 
points are found in diffraction experiments on glasses.^ Quite generally, 
whenever a system is a superposition of two (or more) components such 
that its dynamic quantities, e.g., the dynamic conductivity or a response 
function, have the form described by Eq. (3), isosbestic points are bound 
to occur. This applies in particular to any kind of two- fluid model em- 
ployed, for example, in phenomenological theories of superconductivity and 
superfluidity. There the density of the two components (e.g., the normal 
and superfluid component) depend on temperature while the total density 
is constant: n = n a (T) + re&(T) = const. Properties of the system are then 
described by the superposition of the two components, leading to a special 
dependence of quantities f(T,X) on T and X of the form 

f(T, X) = n a {T)f a (X) + [n- n a (T)]f b (X) . (5) 

This implies the crossing of curves for different temperatures T at a sin- 
gle point X* determined by f a (X*) = fb(X*). Whether the surprisingly 
sharp isosbestic points found, for example, in the optical conductivity of 
Nd2- x Ce x Cu04_j / (Fig. can be explained in this way still has to be 

investigated. 

Disregarding critical phenomena with scaling behavior and the special 
linear dependence given by Eq. (3), isosbestic points in dynamical quanti- 
ties of correlated electron systems have so far only been noticed but never 
explained. As in the case of the crossing of specific heat curves there are 
two separate questions to be answered: (i) Why do curves of frequency- 
dependent quantities cross at all, and (ii) under what circumstances is the 
crossing region confined to a narrow region, or is even point-like? In analogy 
to the entropy sum rule in the case of specific heat curves, a good starting 
point for such an investigation is the study of frequency sum rules. Examples 
are the f-sum rule for the dynamical conductivity and sum rules involving 
the spectral function. 



(4) 
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In this paper we investigate crossing points in frequency by studying 
a particularly basic quantity — the local (momentum-integrated) spectral 
function A(uj, U) of the Hubbard model, a quantity which obeys fundamental 
frequency sum rules. In the following we will write A(u), U) = A(u); the 
parameter U will only be written when explicitly needed. For simplicity we 
will employ the Bethe lattice in the limit of infinite coordination number. 
In this case the self-energy £ becomes local and can be calculated using 
dynamical mean-field theory (DMFT)P E u ' 2 l ' 22 ' 2a ' 24 ' In Sec. EH we discuss 
the behavior of A(co) and in Sec. 12.2.1 the location of crossing points. Their 
behavior at small interaction is studied in Sec. ED We use weak-coupling 
perturbation theory (Sec. l3.lTI) to investigate how the crossing points of A{u) 
depend on the non-interacting bandstructure and the interaction strength U. 
This is performed for both the Bethe lattice (Sec. 13. 2. j) and a model density 
of states (Sec. 13. 3^]) . which differ in their van-Hove singularities. A conclusion 
in Sec. W\ closes the paper. 



2. CROSSING POINTS OF THE SPECTRAL FUNCTION 

2.1. Spectral function in DMFT 

We consider the paramagnetic phase of the Hubbard model © in 
DMFTP | 2U | 21 | 22 | 23 I^I i n this case £ does not have any site- or spin- 
dependence, and the local Green function G{uj) becomes 

G(o;) = Go(a; + M-E(w)). (6) 

The non-interacting Green function Go(z), written as a function of the com- 
plex variable z, is given by the Hilbert transform of the non-interacting 
density of states (DOS) p(e), 

G (z)= [ de^. (7) 
J z-e 

The local spectral function A(u) is then given by 

A(u) = — G"(w + iO). (8) 

7T 

(We use single and double primes to indicate real and imaginary part of a 
complex quantity.) In Fig.[2]the spectral function is plotted for various U in 
the case of half- filling. The numerical data are taken from Ref . 1251 where the 
effective impurity problem was solved with numerical renormalization group 
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co/W 

Fig. 2. Spectral function of the Hubbard model, calculated with DMFT 
for the Bethe DOS [Eq. (jUJ)]; the data are from Ref. 12*51 (a) Sharp crossing 
points (isosbestic points) occur for intermediate U (U /W = 0.8, 1.0, . . . , 2.0). 
(b) For small U the crossing points move up and inward as U is increased 
(U/W = 0.2, 0.4, 0.6). 



(NRG). As usual for model calculations within DMFT the semi-elliptic Bethe 
DOS 

with half-bandwith W was used as non-interacting DOS; we set W = 1 as 
energy scale. 

For all values of U in Fig. [21 the system is assumed to be in the Fermi- 
liquid phase. As the interaction increases, spectral weight is redistributed 
from the vicinity of the Fermi level to the Hubbard subbands which are 
peaked at frequency oj ~ ±U /2, and the well-known three-peak structure 
of the spectral density emerges. The central spectral peak vanishes at the 
metal-insulator transition, which occurs at U = 2.92 VFP^EH 
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2.2. Crossing points 

There are three crossing points visible in Fig. [21 The first one, at uj = 0, 
has a simple explanation: In infinite dimensions Luttinger's theorem en- 
tails that the chemical potential fi in the interacting system must be shifted 
according td^ 

M = Mo + S(0) (10) 

in order to keep the number of particles fixed, where /j,q is the Fermi energy 
for the non-interacting system with the same number of particles. Thus the 
value of the Green function G(uj) = Gq(uj + n — S(cj)) at the Fermi level is 
independent of U in the metallic phase, and the crossing point at uj = is 
exact. That the curves do not really cross but only touch is a consequence 
of the particle-hole symmetry at half-filling, which implies A(uj) = A(—uj). 
For any other density or an asymmetric DOS there is a true crossing point 
at uj = 0. 

The other two crossing points of A(uj) are situated at uj ~ ±1. In 
contrast to the pinning at uj = 0, their explanation is more complicated. 
The crossing region at uj ~ ±1 is very narrow for intermediate values of U 
(0.8 < U/W < 2.0, Fig. Hi). For smaller values of U (Fig. Hp), we observe 
that the region of crossing points is less confined. The crossing points move 
up and inward as U is increased. However, as shown below the crossing 
frequency does not depend strongly on U up to U ~ 2. 

Let us first establish that curves A(u, U) vs. uj for different values of U 
always cross at some frequency uj* , defined by 



8A(uj,U) 



dU 



= 0. (11) 



That two such curves must cross follows from the sum rules 

r0 



71 

dujA(uj,U) = -, (12) 

oo ^ 



dujA(uj,U) = 1, (13) 

where n is the density. The derivatives of these equations with respect to U 
yields 

r0 8A(uj,U) f°° dA(uj,U) 



and thus dA(u, U)/dU vanishes either identically or changes sign at least at 
one point in each of the intervals (— oo,0) and (0, oo). This is the location 
of the crossing point. 
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Fig. 3. (a) The crossing point oj*(U) of the spectral function for the half-filled 
Hubbard model as a function of U and (b) A*(U), the value of A(uj) at this 
point. Each data point corresponds to the intersection of two neighboring 
curves shown in Fig. [21 The dotted lines delimit the range of U in Fig. [2K 
and b, and the shaded bars mark the corresponding width of the crossing 
region Fig. [2^. Dotted curves are the result of second-order perturbation 
theory (Sec. EJ). Due to particle-hole symmetry of the half-filled system it 
is sufficient to consider only u > 0. 



However, the existence of a solution to Eq. (jl 1|) is not sufficient for the 
observation of a sharp crossing point as in Fig. [2^- in general the crossing 
frequency u*, and also the value A* = A(u*(U),U)) of A(u) at this point 
depends on U. This becomes directly evident from Fig.|31 where uj*(U) and 
A*(U) are plotted as a function of U. Since these data were obtained from 
the intersections of curves A(u, U) differing by AU = 0.1 W, while a solution 
to Eq. (jTTj) corresponds to taking AU — ► 0, they have to be understood as 
rather rough estimates for oj* and A*. Furthermore, the resolution of the 
NRG data at co = ±W is already much lower than in the vicinity of the 
Fermi level. Nevertheless, Fig. |31 clearly shows that the strong variation of 
A* with U (for small U) causes the crossing region to become much narrower 
when only intermediate U are plotted as in Fig. |2K- By contrast the value 
of oj* remains close to the non-interacting band-edge uj = W for all U up to 
U « 2W. 
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The variation of A* and w* with U, i.e., the derivatives dw* /dU and 
dA* jdU are not independent. The latter is given by 



dA* dA(w, U) 



dU 



dU 



+ 



dA(w, U) 



dw 



dw* 
~dU 



(15) 



and since the first term vanishes by definition [Eq. we obtain the basic 
relation 



dA* dA(w, U) 



dU 



dw 



dw* 

~~au' 



(16) 



In the present case dw* /dU is small for all U < 2W, while the derivative 
dA/dwlu* diverges as U approaches zero. This is because A(w) approaches 
the non-interacting DOS p(e) [Eq. ©] and w* moves to the band-edge where 
p(e) has a van-Hove singularity. Thus an almost constant crossing frequency 
w* is not in contradiction to a large variation of A*(U). 



3. CROSSING POINTS AT WEAK COUPLING 

As discussed above, well-defined crossing points in A{w) exist not only 
at intermediate U (Fig. |2p.) but also at small U (Fig. (2p) ■ We now investigate 
the behavior of w* in the limit of small U by means of many-body weak- 
coupling perturbation theory (Sec. I3.l!j) . We first show that w* approaches 
the band-edge as w* = 1 + OiU 2 ) if there is a van-Hove singularity at the 
band-edge (i.e., a divergent slope of p(e), see Sec. 13. 2. |) . and later contrast this 
by considering a DOS which vanishes linearly at the band-edge (Sec. 13. 3-1) . 



3.1. Weak-coupling perturbation theory 

We start with a short summary of weak-coupling perturbation theory for 
the local Green function in the paramagnetic phase of the Hubbard model. 
For U <C W one can use a perturbative expansion 

sM = J]f/%H (17) 

n>l 

for the self-energy S, where the functions H n (w) are given by the sum over 
all irreducible Feynman diagrams with n vertices. It is convenient to per- 
form this expansion relative to the Hartree approximationPS The series (j!7|) 
is then rearranged such that each line in the Feynman diag ranis for X] n is 
replaced by the Hartree expression for the Green function, and in turn all 
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Feynman diagrams which contain first-order self-energy insertions are omit- 
ted. Then Si is just given by the static Hartree self-energy Sh = Un/2, and 
only a single Feynman diagram contributes to £2(0;). 

In the following we discuss only the case \x = U/2, for which the Hubbard 
model on a bipartite lattice becomes half-filled and particle-hole symmet- 
ric. This entails a number of important simplifications All odd terms in 
Eq. (|17)l apart from Eh vanish and for the even terms the symmetry relations 
E^(w) = -H' n (-uj) and Tfc(u) = S'^-w) hold, so that E(0) = S H = U/2. 
This is of course consistent with Eq. ifTUl) . 

The calculation of diagrams is considerably simplified in infinite dimen- 
sions, where multidimensional momentum integrals can be reduced to one- 
dimensional integrals over the non-interacting density of states. For further 
reference we mention a convenient expression for the second-order diagram. 
In case of particle-hole symmetry £2(0;) can be written 

/oo 
dXe iXz [2A'(X) 3 - QA"(X) 2 A\X)], (18) 
-00 

where 

/oo 
d£e- iX Zp(Of(0 (19) 
-00 

and /(£) = l/(er^ + 1) is the Fermi function. 



3.2. Bethe density of states 

As a representative of a DOS with a van-Hove singularity at the band- 
edge we first study the Bethe DOS [Eq. ©]. To solve Eq. JTTJ) for uo* we 
insert Eqs. @ and © and obtain 



dA 



dG (z) (d\t dZ(cu*] 



dz V dU dU 



= . (20) 



Here the chemical potential \x has to be differentiated with respect to U 
because the spectrum in (jllj) corresponds to a given density. One can use 
Eq. (jnjj) to rewrite this derivative as dfi/dU = dZ(0)/dU. From Eq. JTUJ) 
and the fact that the first-order self-energy Sh is static, it follows that 
H - E(u>) = 0{U 2 ) for C7 -► 0. Thus dfi/dU - dH/dU vanishes identically. 
In particular one has \i— = — t/ 2 £2(w)+0(£/ 4 ) for half-filling. However, 
instead of looking for a sign change in dA/dU as in Eq. (jllj) one can more 
conveniently investigate the sign changes in dA/d(U 2 ). Eq. (|20j) for the 
crossing frequency at U = 0, 

w = lim lo*(U) , (21) 
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then becomes 
OA 



d(u 2 ) 



\ / Z=OJ* \ / Z=U1* 



For the Bethe DOS the non-interacting Green function Q and its 
derivative dGo/dz are given in analytical form as 

dG n 2z 



G (z) = %z - v^Tv^ + T), -^ = 2--====, (23) 

— ly/z + 1 

where the complex square root denotes the principal branch. When solving 
Eq. (|2Uj) for small U, special care has to be taken to handle the singularity 
of dGo j dz at the band-edges correctly. The leading contribution is given by 



dG 



dz 



for u T 1 

/l— u 
V2 



2 =^+jo [ 77=7 for w I 1 



The perturbative self-energy ^(uj) is continuous and at least differentiable 
once as seen from Eq. (|T5|) and lfTU|) . To obtain cL4/<9(?7 2 ) close to w = 1 we 
can thus replace £2 by its value at the band-edge, £ 2 '(1) = v z ~ —0.2926 
and S' 2 (l) = u 2 ~ -0.0692. The result is 



aA(u;) f Tyfe^ for w|l 
dlP ~ I ^ for w 1 j_ 



(25) 



for C/ — >■ 0. Because V2/112 > this implies a sign change of the derivative 
dA/d(U 2 ) at u = 1, and thus the crossing point tends towards the band-edge 
in the limit of small U. 

Because the self-energy has an expansion only in even powers of U one 
might expect that the same is true also for the crossing point oj*(JJ) as a 
function of U. In particular, this would imply that for small U the crossing 
frequency varies only weakly in the sense that the linear term oj* ~ U is 
absent. On the other hand, the square root in Gq(z) at z = 1 may cause 
to*(U) to become nonanalytic at U = 0. However, as we will now show that 
these nonanalytic contributions occur only in higher order terms, and that 
the solution of (|2Uj) for small U is of the type 

oj* = 1 + aU 2 + Q{\Uf). (26) 

To single out the nonanalyticity of dGo/dz in Eq. 1)20(1 we introduce the 
notation oj* = 1 + 5oj and z = oj* — = 1 + 5z, i.e., Sz = 5oj — E(w*). 

Eq. (J2UJ) then reads 

1= (sz-*)"Ci+ (Sz-vYCn, (27) 
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where we introduced 

Cl = VTTWz du ) I { au ) ' (28) 
° 2 = VTTWz du ) I \ du ) ■ (29) 

Let us now write bz = re^ whence we obtain 

= C2 cos ^ - Ci sin ^. (30) 

With the ansatz (ffi)]) one has y/r = 0(|J7|) and tan</> = bz" jbz' ~ —V2/(a — 
u 2 ) + 0(|Z7|) as U -> 0, whereas Ci = v 2 /«2 + 0(C/ 2 ) and C 2 = 1 + C(C/ 2 ). 
Thus (ESJ) solves Eq. (gOJ) provided that 

u 2 cos — — U2 sin — — ► (31) 

for U — > 0. The solution of this, tan(</)/2) — > V2/U2, together with the 
behavior of tan <fi noted above, yields after some manipulation 



Uo I V. 

-} 1 + - 

I V U: 



1 + 2 = o.653. (32) 
"4/ 

The results of this section are easily generalized for any non-interacting 
DOS with a van-Hove singularity p(oS) ~ — W\ at the band-edge u = 
dbW, as, e.g., for the DOS for a simple cubic lattice. This is because the 
dependence of Gq(z) = J dep(e)/(z — e) and dGo/dz = — J dep(e)/(z — e) 2 
on p(e) is linear. The singular part of dGo(z)/dz can be split off and the 
remaining continuous part does not enter the first-order result for cu* . We 
thus conclude that a van-Hove singularity at a band-edge implies that in its 
vicinity a crossing point in the spectrum exists for small U, provided that 
£2(1) an d ^2(1) nave the same sign. 



3.3. Model density of states 

The situation is different if there is no van-Hove singularity at the 
band-edge of the non-interacting DOS, e.g., if the DOS vanishes linearly. 
In this case its derivative at the band-edge does not diverge. Then the 
weak-coupling crossing point loq need not be located at the band-edge. We 
now demonstrate this by investigating the family of model functions 

Pa(e) = P^(l-\e/Wn&(W-\e\), W = I, a = 1,2,3,... (33) 
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as non-interacting DOS. Here 0(x) is the step function and again the half 
band-width W = 1 has been chosen as energy unit. 

The functions p a (e) are all linear at the band-edges but become increas- 
ingly steep as a increases (Fig. |3Ji). In the limit a — > oo p a (oj) approaches 
the box-shaped DOS Poo{u) = 0(1 — |u;|)/2, for which a similar calculation as 
in the last section shows that u>o = ±1. We therefore write uj^ = 1 — A a /a 
for the crossing frequency in the limit U — * for the DOS p a {e). As shown in 
the following, A a — > Ao for a — > oo, i.e., there are corrections of the crossing 
frequency with respect to the band-edge, but these corrections become small 
as the DOS develops a discontinuity there. 

We first write Eq. (J22J) as 



(dGo/dz)' 



(dGo/dz)" 



z=l-A a /a+iO (E^) 



(34) 

z=l-A a /a+iO 



The second-order contribution to the self-energy is shown in Fig.0Jo. Similar 
to the last section, the nonanalytic behavior of dGo jdz at z = 1 determines 
the behavior of the solutions of this equation. For A a <C a and a —> oo the 
leading contribution to the derivative [see Eq. (J7J)] 

dGa a + l fl dx 1 = X ° (35) 



dz 2a J_ x (1- A« +i0 -x) 2 ' 



is given by 

cfGo a f 1 



dz 2 \A a J (x-A a + iO) 2 



+ Re / dx- — j + ivre" Aa . (36) 
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Fig. 5. Location of crossing points in the spectral function in the limit 
U — > for the non- interacting DOS p a (e) [Eq. (f33|> ] . The triangles mark 
the numerical solution of Eq. the dotted lines shows the lowest order 

approximation [Eq. I|37|)]. 



On the other hand the right-hand side of (|34|) is finite for a —* oo. We can 
thus analyze the behavior of the solutions for large a by letting A a — ► and 



a 



oo in E 



(a) I 



|l-A«/a- 



To lowest order we then have 



(a) 



1 - A /a, 



(37) 



where Aq is the a solution of 



A f 



+ Re 



dx- 



A + iO) 2 



7T 



( E (oo)y/ 



(SH)' 



U)=l 



(38) 



Since S^ 00 ^ (w) can be calculated directly from Eq. (JT5|) and (fTU]) we can solve 
(jSHj) numerically. The result, Ao = 0.0750, is compared in Fig. [5] to a direct 
numerical solution of of Eq. Q18j). and which is possible for not 

too large a. Fig. |S] clearly demonstrates that in the limit U — * the crossing 
points are not located at the band-edges if there is no van-Hove singularity. 
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4. CONCLUSION 

In this paper we presented an analytical investigation of isosbestic points 
in the local spectral function A(u, U) of correlated electrons described by 
the Hubbard model. The existence of such crossing points in the curves of 
A(lo, U) vs. uj for different values of the interaction U is due to sum rules for 
A(u,U). For small U the frequency corresponding to the isosbestic point 
may be derived within weak-coupling theory. In particular we showed that 
if the non-interacting DOS has a van-Hove singularity at the band-edge, 
there is in general a crossing point at uj* = W + 0(U 2 /W) for U — ► 0. By 
contrast, if the DOS vanishes linearly at the band-edge the crossing frequency 
is located at u*(U = 0) ^ W. 

We established that, in general, the crossing frequency uj*, and thus the 
value of the local spectral function at the crossing point, A(u)*), sensitively 
depends on the form of the non-interacting DOS p{e) at the band-edges. This 
is quite different from the near-universality of the crossing point values of the 
specific heat at high temperatures.^ We therefore conclude that the study 
of isosbectic points in dynamical quantities can provide valuable information 
about microscopic energy scales of strongly correlated electron systems. 
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